Random tilings of high symmetry: I. Mean-field theory 
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We study random tiling models in the limit of high rotational symmetry. In this limit a mean-field 
theory yields reasonable predictions for the configurational entropy of free boundary rhombus tilings 
in two dimensions. We base our mean-field theory on an iterative tiling construction inspired by the 
work of de Bruijn. In addition to the entropy, we consider correlation functions, phason elasticity 
[ and the thermodynamic limit. Tilings of dimension other than two are considered briefly. 

> ' 

^ ' I. INTRODUCTION 

a. 

The discovery of quasicrystals by Shechtman, et al^, and the proposal of Levine and Steinhardt- of their possible 
relationship with Penrose tilings^, motivated widespread investigation of tilings as models for quasicrystal structure. 
It is presumed that energetically favored atomic motifs form geometrical tiles, and that these tiles may be arranged 
quasiperiodically in space to describe the quasicrystal structure. We consider tiling models defined as projections 
!— 1 1 from a higher D-dimensional lattice into a lower d-dimensional physical space. For example, 10-fold symmetric tilings 
may be projected from D = 5 into d = 2, and icosahedral tilings may be projected from D = 6 into d = 3. Figure^ 
displays examples of tilings constructed as projections from D-dimensional simple cubic lattices into d = 2-dimensional 
i physical space. The difference between the higher dimension D and the tiling dimension d is called the codimension. 
' An important unsolved problem is understanding the interactions among tiles that may favor particular tile ar- 
l/"") , rangements over others. The suggestion by Elser 4 that random tilings spontaneously exhibit quasiperiodicity, created 
a subfield within the area of tiling theorj^. It was showri^, for at least one atomistic model quasicrystal^, that 
quasicrystalline order emerges with random tiling, rather than Penrose-like, order. The best description for real 
quasicrystalline materials remains an open problem, depending whether one assumes that energetic interactions favor 
| quasiperiodicity or that entropy stabilizes the quasicrystal structure at high temperature o 5 ' 6 ' 9 . 

In addition to their role in the theory of quasicrystals, random tiling models appear in a number of other interesting 

S scientific contexts. Their combinatorial properties, for example their relationship to generalized partitionsiS*iiii* make 
them interesting models for study within pure mathematicsi2ii2ii£ii£iiLi2ii2i2S. Some random tiling models provide 
examples of exactly solvable models2ii22iS of interest within statistical mechanics. We note a relationship with 
algorithms for sorting Iists2i2^. Random tiling models have also been proposed as models for elastic membranesii*2&. 

The random tiling theory of quasicrystals focuses on two important properties: The tiling entropy contributes to the 
configurational entropy of the quasicrystal, reducing its free energy and enhancing thermodynamic stability against 
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FIG. 1: Examples of free boundary d — 2 tilings with D = 3, 5, 7 and 9. 
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TABLE I: d = 2 strain-free rhombus tiling entropies. 



D 


2 


3 


4 


5 


oo 


CT 





0.323 


0.434 


0.481 


0.568 



other competing phases 9 ; Variation of the entropy with average phason strain defines the phason elastic constants. 
Table [I] summarizes existing data on the configurational entropies of d — 2 rhombus tilings, showing a trend of 
monotonic increase with D. Table ITT1 shows corresponding data for d = 3 rhombohedron tilings. 

TABLE II: d = 3 strain-free rhombohedron tiling entropies. 



D 


3 


4 


5 


6 


oo 


CT 





0.214 


NA 


0.24 


< 0.323 



Techniques employed in theoretical studies include analytic exact solutionis, Monte Carlo computer simulation^! 
and numerical transfer matrix methods^. While these methods prove successful at calculating quantitative values 
for entropy and elastic constants, they are not always easy to implement, and they often do not enhance our intuitive 
and qualitative understanding of the problem. The technique of mean-field theory^ in statistical mechanics provides 
an approximation that can be analytically simple to implement and can provide a direct link between quantitative 
calculation and our qualitative picture of the physics of the problem. Mean-field theory often becomes exact in certain 
limits, such as the limit of infinite dimension. A mean-field theory for random tilings has never been precisely defined. 

The search for an analogue of mean-field theory applicable to random tilings motivates our study of random tiling 
models with high rotational symmetry. We anticipate^ that the infinite D limit may prove easier to analyze than 
specific finite values of D. Experience gained in the study of this limit may later be transfered back to the finite D 
values of greatest physical interest. Furthermore, a study of D dependence of entropy or elastic constants might in itself 
yield qualitative understanding of the values of entropy and elastic constants for small D. Other researchers^SiiS 
have studied deterministic quasiperiodic tilings of arbitrary rotational symmetry, but have considered neither the high 
symmetry limit nor the case of random tilings. 

We mention, in addition, that high symmetry tilings are good candidates of models for topologically disordered 
graphs, different from the usual random graphs^ because the coordination number of a vertex is a random variable, 
with algebraically decreasing long-range correlations. Quantum diffusion on such tilings has recently been studied by 
Vidal et al. and was found to exhibit an interesting diffusive behavior™. These tilings could even prove interesting 
as models for structural glasses. 

A word of caution should be given on the title referring to "high symmetry" tilings. In two dimensions, this 
statement is unambiguous owing to the possibility of generating tilings whose Fourier spectrum will display point 
symmetries of arbitrary order. In three dimensions however, the highest point group symmetry, the icosahedral one 
of order 120, is already reached with codimcnsion 3 tilings. "High symmetry" , in three and higher dimensions, has 
therefore a less formal meaning in our description. For a given tiling in dimension d, we can map the whole set of 
edge orientations onto points on a (d— l)-dimensional hypersphere. As D increases, these representative points cover 
more and more densely the hypersphere, approaching, but never meeting, the highest continuous point symmetry of 
the hypersphere itself. This is what is meant by "high symmetry" in the title. 

Two previous papersS^*^ presented initial studies of this problem. The first paper proposed an upper bound on 
the entropy in the limit of large D and discussed problems associated with the thermodynamic limit of tiling models. 
The second one^ presented a preliminary mean-field approach of the entropy calculation. In the present paper we 
explain the previous work in greater detail. In addition, we inspect spatial correlations and phason elasticity. All 
tilings considered in this paper are free boundary tilings (see section [HJ . Fixed boundary tilings are analyzed by 
computer simulation in the following paper—. Table ILTTI summarizes the available estimates of the large D entropy. 

The organization of this paper is as follows: We begin in section^with definitions and concepts that apply generally 
to random tiling models. Then, in section II I II we describe our mean-field theory of two-dimensional tilings of high 
rotational symmetry. We start with an upper bound fsection |lII If) then a more accurate estimate fsection |lII 21 based 
on path counting arguments. We investigate the nature of spatial correlations in section llll 3l followed by consideration 
on finite D corrections. Section I1VI tackles de Bruijn line interactions and phason elasticity. Tilings of dimension 
d 7^ 2 are examined in section Ivl 
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TABLE III: Entropy estimates in the large D limit. 



Value 


Method 


Reference 


0.998 


Rigorous upper bound 




0.693 


Conjectured upper bound upper bound 


eq. JTTTTft 


0.6 


Extrapolated exact values 


Ref. f2&] 


0.598 


Mean-field theory 


eq. HIII.14fl 


0.568 


Monte Carlo simulation 


Ref. p£] 


0.231 


Rigorous lower bound 





II. CHARACTERIZATION OF TILINGS 

In this section, we introduce the basic material used throughout the paper to characterize and manipulate two- 
dimensional random tilings. 

A tiling is a covering, without gaps or overlaps, of a given region of a d-dimensional Euclidean space. In the present 
paper, the tiles are d-dimcnsional rhombohedra, which we will generically call "rhombi" in the following. 

Tiling systems can have different boundary conditions, such as free, fixed or periodic ones. For example, the tilings 
in figure ^ illustrate free boundary tilings of various rotational symmetries. The free boundary thermodynamic limit 
is taken as follows: Consider all tilings that cover a circular region of space 1Z (whose area will tend towards infinity) . 
For a given TZ, tilings contain fixed numbers of tiles of each type, with relative frequencies governed by the average 
composition of the tiling we wish to study. The total number of tiles is chosen so they encompass a total area slightly 
larger than area of the region TZ that must be covered. Every tile must intersect TZ or share an edge with a tile that 
does so. This ensures a reasonably compact tiling and gives us needed control over the composition. 

It is generally supposed that free and periodic boundary tilings have the same entropies (provided their compositions 
match) at the thermodynamic limit. But fixed boundary tilings display a less usual behavior: such boundaries have 
such a macroscopic effect on tilings that their entropy per tile at the infinite size limit can be strictly smaller than 
the free boundary one. For example, in the case of "hexagonal" tilings of 60° rhombi (D = 3 and d = 2), when all 
orientations of rhombi occur with equal probability, the fixed boundary entropy per tile a fixed — 0.261 (— ) while the 
free boundary entropy per tile cr/ ree = 0.323 ( 3S ). In addition to a qualitative argument by Elser— , this phenomenon 
has more recently been rigorously explained and describecUii&Si: the local entropy density can be calculated in any 
point of the tiling; it displays a gradient between the regions near the boundary, where the entropy is vanishing, and 
the central region, where the entropy is a free boundary one. In other words, it is only at the very center of the tiling 
that the tiling loses the influence of the boundary. 

1. Membrane representation of tilings 

Here, we briefly recall how tilings are conveniently coded by (mono- valued) functions from the "physical real space" 
to the "perpendicular" space. 

Rhombic tiles can be considered as the projections onto a d-dimensional space of the d-dimensional faces of a 
hypercube in a space of higher dimension D. Conversely, any rhombic tiling in the d-dimensional space can be "lifted" 
to a d-dimensional continuous membrane embedded in the Z D hypercubic lattice. When this membrane is projected 
back into <i-dimensional space, the facets of the Z?-dimensional lattice project precisely onto the d-dimensional rhombic 
tiles. Such a membrane is said to be continuous and directed to emphasize the fact that, when projected, is does 
not exhibit any gaps or overlaps. We say that we are dealing with D — > d membranes or tilings. This membrane 
representation of rhombus tiling 

,5.11.26.39 ig a 

useful tool in quasicrystal science. 
Notation related to the D-dimensional geometry and projection is described in appendix 1X1 Let (e^^i^..^ denote 

the basis of the Z D hypercubic lattice and denote by e| the normalized projection of §j on the (i-dimensional space 

where tilings lie. The projection direction is (arbitrarily) chosen so that the basis vectors ef all have equal length and 
angles that are multiples oiir/D (see appendix lAl. We call this collection of vectors a regular fan. 

A rhombic tile has d edge orientations {ef , e|, • • • , e|} with (1 < i < j < ■ ■ ■ < k < D). We denote this tile species 

as Tij...k- Note that there are different species of tiles. Examples of D — > 2 tilings are displayed in figure^ The 
first one can be seen as a membrane embedded in a 3-dimensional cubic lattice, as may be more clear in the following 



figure [21 
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FIG. 2: Three-dimensional image effect. Gray rhombi can be seen as horizontal terraces. 

Such random tilings have internal degrees of freedom, the so-called localized phason flips, which consist of local 
rearrangements of tiles. Such local flips are displayed in figure |21 In d dimensions, a flip involves d + 1 tiles that fill 
a small zonotope^. Note that in 2 dimensions, the ergodicity of tiling ensembles via elementary flips is proveniSii4, 
whereas it is an open question in 3 and more dimensions, even if some significant progress has been achieved recently^. 
This local degree of freedom is the source of configurational entropy in random tilings. 




FIG. 3: Two examples of elementary flips in 2 and 3 dimensions, involving respectively 3 and 4 rhombic tiles. 

A functional representation of these directed membranes will prove useful. Denote the ci-dimensional "physical" 
or "real" space containing the actual tiles by £" = R d , and the perpendicular (D — d)-dimensional space by S 1 - = 
R D ~ d . Since membranes representing tilings are directed, they are the graphs of single-valued continuous functions 
h : £ " — > £-*-. Such functions h are called the height functions of the corresponding tilings. The detailed construction 
of the spaces £" and £ ± is described in appendix lAl when d = 2. Executing an elementary flip of vertex v changes the 
value of the height function locally at vertex v. 

Since the membranes are embedded in a hypercubic lattice, they are irregular and faceted at on the length scale 
of individual tiles. The local height function h is usually coarse-grained to get a smoother function h which closely 
tracks the local function h over long distances but which is nearly insensitive to individual single vertex flips 5 ' 11 - 42 . 
The function h describes a smoothed membrane. Then one defines the phason gradient E = Vft. This (D — d) x d 
dimensional tensor controls the fractions of the different species of tiles. 

The entropy can be written as a functional of the phason gradient 5 . Moreover, the random tiling model hypothesis 
states that this latter entropy has a unique maximum, corresponding to tile fractions maximizing the symmetry. The 
orientation of the real space £ is chosen so that the gradient is zero at this maximum, and the model states that the 
entropy density has a quadratic behavior near this maximum: 

(J = cr max — i Vh ■ K • Vh + o(\Vh\ 2 ). (Ill) 

By analogy with an usual clastic theory, the tensor K is called the tensor of phason elastic constants. 

The best known rhombus tilings of the plane is the Penrose^ tiling, which displays an exact quasiperiodicity. The 
membrane representation of a Penrose tiling is a flat, horizontal, sheet. The phason strain of this tiling vanishes. 
In addition, the Penrose tiling is virtually unique (up to translations) because it forbids localized phason flips. The 
random tiling hypothesis implies that the most typical tilings in a random tiling ensemble mimic the long-range order 
of the Penrose tiling, while enjoying the entropy of short-range fluctuations. 
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Random and perfect tilings with vanishing phason strain exhibit the highest symmetry consistent with the projection 
from D — > d. In two dimensions, this symmetry includes^? 2_D-fold rotational symmetry. A given tiling itself generally 
has no point about which it is invariant under rotations, but two-point and higher correlation functions do exhibit 
the symmetry. In addition, there are two families of mirror symmetries, denoted "along" and "between" according to 

the relation between the mirror plane and the basis vectors el . 

2. De Bruijn dualization and iterative construction of tilings 

The de Bruijn dualization 4 ^ provides another representation of random tilings. We present this technique in two 
dimensions and generalize it to higher dimensions when needed. 

De Bruijn g^ida^^^ are made up of lines, the de Bruijn lines, that are sometimes called worms. These lines are 
made of adjacent tiles sharing an edge of given orientation. It is always possible to extend these lines through the 
whole tiling up to a boundary tile. An example is presented in figure 0] Since there are D different edge orientations, 
there are D different de Bruijn families. A rhombic tile corresponds to the intersection of two lines of different families. 

Hence there are = D(D — l)/2 tile species. There are no triple intersections, and lines of a same family never 

intersect, even in an infinite tiling, because no rhombus has four parallel edges. 




FIG. 4: Iterative process for the construction of D — ■> 2 tilings. Left shows paths on a re-oriented 3-^2 tiling. Right shows 
the resulting 4 — > 2 tiling. A non-crossing contact between two de Bruijn lines of the same family is emphasized. 

This de Bruijn representation of tilings suggests an iterative construction of D — > 2 tilings in terms of directed 
random walkers on D — 1 — > 2 tilings. The basic idea is illustrated in figure 0] A D — I — > 2 tiling is sheared to make 
room for a new de Bruijn line family by taking the vectors e| that initially lie at angles iri/(D — 1) and replacing 
them with angles ni/D (see Fig. left). Next, paths are chosen on the D — I — > 2 tiling (dark lines in Fig. left). 

They go from bottom to top, as symbolized by a "time arrow" t. These paths must be "opened" in the direction e^ D 
in order to form de Bruijn lines (shaded in figure 0] right) of the D-th family and therefore a D — > 2 tiling. These 
paths are directed in the positive time direction and they do not intersect within the same family. The figure shows 
a contact among paths that is not an intersection. Conversely, de Bruijn lines of family D in a D — > 2 tiling may be 
collapsed to directed walks in a D — I — > 2 tiling. There is thus a one-to-one correspondence between D — > 2 tilings 
with p de Bruijn lines of family D and collections of p non-intersecting walks on D — 1 — > 2 tilings. 

3. Grid parameters, phason strain and tile fractions 

In this subsection, we introduce several macroscopic quantities that characterize a two-dimensional random tiling. 
Apart from the phason gradient E previously defined, a tiling can also be characterized at large scale either by the 
tile fractions or by geometric parameters. This subsection also provides the relationships between these quantities 
(|11 3 20 as well as useful relations that constrain them l|ll 3 31) . 
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1. Definitions 

We now address two macroscopic characterizations of de Bruijn lines, the mean line spacings and orientations, and 
link these to the phason strain E and the tile fractions n« (the fraction of tiles that are of type T^). De Bruijn lines 
of family i are characterized by the mean spacing Zj between two lines of the family i and the mean angle fa that 
lines of this family make with their nominal direction . The nominal direction (denoted by Si) is normal to ef (see 
figure EJ). For a maximally symmetric tiling with vanishing phason strain E = 0, rotational symmetry guarantees 
that the mean line spacing U is independent of the line family i. We shall calculate its common value shortly. The 
strain-free tiling's mirror symmetries guarantee that each angle fa = 0. Deviations of U and fa from their strain-free 
values control nonzero values of E and determine the tile fractions n, 3 . By definition, li > and — 7r/2 < fa < tt/2. 

Figure |S] represents lines of the i-th family as sequences of tiles Ty , each tile defined by the vector ef and some 
other vector ejj . Let y be the relative position of a given tile in a de Bruijn line, and follow this line in the plane. 
The position y m — y of the m-th tile is the sum of the vectors defined by the successive edges of the tiles TJy of the 

line between tile and tile to. These vectors are a priori equal to iejj. Therefore we define f, = ±ejj where the sign 
is chosen so that fj • Si > 0. For large to, the number of tiles T^- on any typical line is proportional to the tile fraction 
riij . The mean direction of this line is therefore 

D 

3=1 

In general, a random tiling ensemble is completely characterized by {h} and {fa}. Since D — 2 families of lines must 
be added to an original square lattice (D = 2) to define a D — > 2 tiling we need 2(D — 2) independent macroscopic 
parameters to characterize a tiling: li and fa for i = 3, 4, ...,D. Note that 2(D — 2) is precisely the number of 
coefficients of the global phason gradient E = Vh. By comparison, in reference faS], 4^2 tilings are characterized 
by 6 tile concentrations dij, and 2 exact relations constrain them. Therefore these tilings are characterized by 4 
independent parameters, as predicted above. 




FIG. 5: Definition of quantities U - the mean distance between two neighbor lines of the i-th de Bruijn family - and fa - the 
mean angle that the lines of this family make with Si. The vector Uj is given by the vector ef rotated by the angle fa. The 
orientation di of a line of the family i is given by the average of the vectors f, which define its tiles. 



2. Relations between macroscopic parameters 

Consider now the tile fractions riij and their relation to {li} and {fa}. The tile number fractions are proportional 
to their density dij (number per unit area) via 

»■-- dl3 (II.3) 



J2k<l d kl 

Conversely, 



n ; 

J2k<i n ki sin|7r(Z - k)/D\ 



dil = ^ Z ±hn i^THi ( IL4 ) 
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because the area of a tile equals sin \9*j | where 6*j — ir(j — i) / D is the oriented angle between e" and . 

Note that 0*- is also the angle between de Bruijn line directions i and j in a strain-free tiling. Generally speaking, 
throughout the text, starred quantities will refer to unstrained tilings with E = 0. 

Relation (|II.3(I is useful because the tile densities are given in a simple fashion by 42 - 46 

djj = \rrii x m.j\, (H-5) 

where m,i (see appendixEJ is the parallel space gradient of the i th component of the smoothed membrane hyperspace 
coordinates. We also define a vector Ui for each line family: it is the unit vector along the shortest direction from 
one line of the family to the next; it is perpendicular to di. Examining figure if one travels a distance U in the 
direction Ui, then one crosses one de Bruijn line of the family i on average and the i th component of the hyperspace 
coordinates increases by one on average. In contrast, if one travels any distance in the orthogonal direction d iy then 
the i th component of the hyperspace coordinates remains constant. Thus 

mi = yUi (II. 6) 



and equation (|II.5|I becomes 

d ij = JT \ u i x u j \ = JT I sin % I > ( IL7 ) 
where the angle between de Bruijn line directions 

On = etj + & - ( IL8 ) 

Next we relate the distances {k} and angles {<fri} to the phason strain E. Owing to relation (|A.10|I of appendix 1X1 
we have 

D 1 

E = ^2~ (g>u z . (II.9) 
Defining E 2 as the sum of the squares of the components of E, then 



E2 = EE 777^ (11.10) 

i— 1 j — 1 J 

We prove in appendix 1X1 that if i ^ j, then ef- ■ ej- = —2cos8*j/(D — 2). Hence 

2 



E 2 



£ _L cOS%COS0*,+£ 4 (11.11) 



D - 2 ^ /,/, ' « ^ I 



To determine the tile fractions n*^ and the mean line spacing I* for strain-free tilings (we recall that starred quantities 
refer to unstrained tilings), remind that the tiling symmetry guarantees U is independent of i, and (pi — 0. Inspecting 
eqs. (|II.7ll and l|II.3|) . we may write riij = C\ sinf^* |. To determine the normalization constant C, note that 

i = 5>£ = c Yl ( n - 12 ) 

i<j 0<i<j<D 

In the limit of large D one then finds C ~ tt/D 2 , so that 

4 = -^|sin0*.|. (11.13) 
The fraction rij of tiles belonging to a given family i of lines, in a maximally symmetric tiling at large D, 

n a - 7^7 = n (IL14) 
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allows us to compute the mean distance between de Bruijn lines of a given family. Since tiles of family i have, on 
average, the same area as tiles of the whole tiling, the above numerical tile fraction is also an area tile fraction. Since 
a fraction 2/D of tiles belong to lines of a family and the width of such a line is 1, the mean distance between lines is 

I* = D/2. (11.15) 

Setting k = I* and Oij — 8*j for all i and j in we find again in a consistent way E 2 = 0. Moreover, eq. (Ill . T|) 

implies 

4 4 

d ij = Jyi I sin °*i I = ~ n ij ( IL 16 ) 

in the large D limit. 

3. Miscellaneous constraints 

Geometrical considerations constrain the parameters {7,} and {&}. We describe a few such constraints here. The 
condition 

4>i-i ~^<<Pi< <Pi+i + J (H-17) 

imposes a certain regularity to the angles <f>i as a function of i. This constraint holds because lines of two consecutive 
families i and i+ 1 intersect only at a tile Tj^+i, with edges ef and e|' +1 as displayed in figure [B] The oriented 
angles (djdj+i) and (<$j<5j+i) have the same sign. Indeed, were these angles of opposite sign, no tile type could match 
at their intersection. As a consequence (djdj+i) = (5i5i+i) — <f>i + 4>i+i > 0. Owing to (8i8i+i) — tt/D, one gets 
4>i+i — 4>i ^ —tt/D and condition l|II.17|) (see figure EJ. Such conditions mean that at large D, when one adds a 
new family of lines in the iterative process, the angle 0_d is highly constrained by the pre-existing underlying tiling. 
For example, if a family i = D is added in a strain-free tiling where <pi = 4>d-\ = 0, the condition (|11.17|l reads 
-n/D < <t>D< n/D. 




FIG. 6: Left: Intersection of two lines of two consecutive families i and i + 1. Right: The corresponding directions Si, di+i 
(outer thick lines) and di, dU+i (inner thin lines). The vectors di and dj+i are constrained by the condition (didi+i) > 0, 
which constrains in turn the angles cj>i and 4>i+i. 

Since amongst the 2D parameters {li} and {4>i}, only 2(D — 2) arc independent, there exist 4 relations between 
them. Of these four we will need this identity 

]T^|sin%||sin^| = l (11.18) 

which reflects the fact that the total contribution of tile areas per unit area is by definition equal to 1. In other words, 
J2i<j \ sm ^i*jl — 1 where is given by eq. I|II.7|) . This relation can also be derived from eq. (|II.4|) . 

In the large D limit we have infinitely many tile species, so we expect that all — > 0. When taking the large D 
limit we wish to ensure that no finite fraction of tile species dominates, with tile fractions that are large compared 
to the remainder. Indeed, such a tiling would essentially be a finite D tiling, with a negligible number of defect tiles 
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added in. Thus we introduce the notion of "bounded" relative fractions of tiles, which proves useful in discussing 
large D phason elastic constants. Consider tilings in which all tile fractions tend to as D — > oo, all vanishing at the 
same rate. Specifically, assume there exists a > such that 

n^^an*^ (11.19) 

holds uniformly for all i and j , which ensures that all tile fractions have the same order of magnitude as in strain-free 
tilings. Owing to eq. 1)11. 4|) . we could alternatively assume that there exists b > such that 

<!,, > b,i; r (11.20) 

Indeed, X)fc<; n ki — 1 thus X)fc<2 n ki sm l 71 "^ — k)/D\ < 1 and dij > riij. Furthermore riij > a n*^ — -kcl/A d*j and 

b = Tra/4. 

In addition, conditions (|II.19(I or 1)11. 20|) constrain the distances h and the angles Oij so that 

k<c I* (11.21) 

where c > is a finite constant, and 

| sin % | > d |sin0?.| (11.22) 

with d > 0. Indeed, these two conditions together with eq. l|II.7fl imply d^ > (d/c 2 )d*y 

Comparable constraints can be stated in the membrane viewpoint, in which case the constraints appear as extreme 
allowed values of E. Using (|II.llfl and constraint l|II.21|) . one finds that the phason strain is vanishingly small at large 
D, 

^ 9 Const. , 

E < — - — . (11.23) 

The only way to achieve finite strain at large D is to allow a finite subset of the tile species to dominate. This large 
D behavior is related to our choice of normalizations of basis vectors both in and in £K Our choice is consistent 
with current practice for finite D tilings. 



III. MEAN-FIELD THEORY, CORRELATIONS AND FINITE D CORRECTIONS 

This section applies the iterative construction of section ITT~2l to develop a mean-field theory of strain-free random 
tilings. We first propose a close upper-bound of log 2 for the entropy per tile, which can be seen as the order of 
the mean-field theory developed below. We then refine the estimate by calculating, within mean-field theory, the 
distribution of tile vertex types. Finally we explore spatial correlations and related issues. 



1. First entropy estimate 



Here we set the upper-bound of log 2 for the entropy. To begin with, we mention that rigorous bounds of the large 
D entropy are known in the case of fixed-boundary tdings of unitary side length (see paper II2&, section II): 

- log 2 < o£ xed < 1.44 log 2. (IILl) 
3 

Since it is demonstrated in the same reference that this entropy per tile is equal to the free-boundary large D entropy 
<Joo, the same bounds should hold for a^. Our first mean-field estimate log 2 belongs to this interval, as well as the 
refined mean-field estimate of the next section. 

Let Pd(Nd) denote the number of ways to follow a directed walk of family D on a D — > 2 tiling. is the total 
length of a single walk, or the cumulative length of a few walks. We will neglect contacts among walks of family 
D, which is justified since i, ~ D. To calculate a typical value of Pd(Nd), consider the problem of constructing 
bottom-to-top paths on a typical tiling. At each step, the path may follow one or more routes (Fig. 0). The only 
requirement is that the path segment leading out of a vertex must contain an upwards component. We estimate the 
number of paths as the product over vertices v of the number of choices N c (v) to be made at each step 



Pd(N d )^Y[n c (v). 

v=l 



(HL2) 
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FIG. 7: A random path (thick line) on a tiling as a succession of vertices with choices (dashed edge = choice not taken). There 
are 5 choices at vertex A (i.e. N C (A) = 5), 2 choices at vertex B and only 1 choice at vertex C. 



Evaluation of equation (|III.2I) requires the distribution of values N c (v) along paths. Because we do not know this 
distribution (yet), we settle for an estimate that yields an upper bound on Po(A r £i), and eventually on the entropy per 
tile Sd- Note that the product in equation (|III.2|I is the Nu-th power of the geometric mean of N c (v). The geometric 
mean of any set of positive numbers is bounded above by the arithmetic mean, reaching this bound only when all 
values are equal. The arithmetic mean N c = Z/2, with Z the mean coordination number (or degree), because on 
average half the tile edges at each vertex have a vertical component in the northerly direction. From Euler's theorem 
applied to infinite rhombus tilings, we know Z — 4, so that N c = 2. We deduce the D-independent upper bound 

Pd{Nd) <2 Nd - (HI.3) 

Let Bd be the number of D — > 2 tilings (free boundary tilings of some size TV). Since each D + 1 — > 2 tiling is in 
one-to-one correspondence with random walks on a D — > 2 tiling, 

B D+1 ^P D (N D )B D . (III.4) 

We iterate this relationship to write 

D-l 



Bd = IT P D'(N D ,). (III.5) 



D'=3 



Inserting relation (|III.3fl for the path counting, we find 

B D <2 N3+ - +Nd -K (III.6) 
Since the number of tiles is N = Y2d' Nd', we obtain the entropy per tile 

<j D = < log 2 = 0.693. (III.7) 

The limiting entropy (Too is bounded above by log 2. 



2. Refined mean-field calculations 



The entropy value log 2 just derived is not exact for two reasons. First, to properly compute the geometric mean 
of N c we need the full vertex path choice probability distribution tt(N c ). Second, our estimate eq. HIII.2|) neglected 
correlations among the numbers of choices at different steps along the path. In this subsection we address the first 
point. In the next one we shall address briefly the second point. 

Following the logic in section ITlI II the limiting entropy is 

aoD = Hm Hm l ^fER , (Hl.8) 

D—*oo k — *oo k 
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where Po(fc) is the number of k-step paths on a D — > 2 tiling. The mean-field approximation assumes that the steps 
of such paths are uncorrelated so that the number of paths is given by eq. i|111.2fl . The vertices v belong to a random 
D — > 2 tiling. Therefore the numbers of choices N c (y) are distributed according to a probability distribution ttd(N c ). 
When D tends to infinity this distribution tends toward a limiting distribution, denoted by 7r(JV c ). Thus, after a 
short calculation, equations l|III.8(l and (|III.2I) become: 

oo 

*% F = E <N c )]ogN c . (III.9) 

N c = l 




FIG. 8: A path going through a vertex v of a D — » 2 tiling gives birth to two vertices ui and V2 of a D + 1 — * 2 tiling. The 
path follows any leg of v with probability 1/g and any arm with probability 1/p. 

To get the mean-field distribution ir(N c ), let us first denote by 7ro(q,p) the fraction of vertices on a D — > 2 random 
tiling with q in-coming edges ("legs") and p out-going ones ("arms") (see figure |SJ). Neglecting correlations, a (q,p)- 
vertex v of a D — > 2 tiling will be visited by paths with probability iTD{q,p) and will then give birth to two vertices 
Ui and i>2, as illustrated in figure |H1 These latter vertices belong to a D + 1 — > 2 tiling. Each leg and each arm of u 
will be chosen with probabilities l/q and 1/p, respectively. The new edge from v\ to vi is arbitrarily chosen to be 
oriented upward in the D + 1 — * 2 tiling and is therefore an arm for v\ and a leg for V2 ■ 

Indeed, a vertex with p' arms and q' legs can be either the "left-son" v\ or the "right-son" vi of a "father-vertex" 
v in a D — > 2 tiling (figure 0. The two son vertices together inherit all their father's arms and legs, plus an extra 
arm and leg (for the path along which the vertex is split). Finally, in addition, the left son gets an extra arm and the 
right son an extra leg from the new bond that joins them. Consequently the left-son has more than one arm and the 
right-son has more than one leg. 

As a consequence, in this mean-field approximation, the probabilities ttd+i(q>p) can be written as linear combi- 
nations of the probabilities t:d(q,p)- The corresponding linear operator will be denoted by A : -kd l— * Td+i- It is 
infinite-dimensional since q and p can be arbitrarily large when D goes to infinity. The operator A is defined by: 

i i i\ ft x \ ^oiq.p) , x , v-» n D (q,p) mum 

TT D+1 (q,p) = (l-5 q , A ) ^ 2gp + (1 - E 2gP " ( ^ 

q>q' — 1 ' — 1 

The first and second term refer, respectively, to left- and right-sons (vi and t^)- The factors of 1/2 expresses that a 
vertex is either a left-son or a right-son with probability 1/2. Moreover, each arm-leg pair is chosen with probability 
1/pq among all possibilities (see figure |SJ). The limiting distribution Tr(q,p) is the fixed point of A, the eigenvector 
with eigenvalue 1. 

The diagonalization of A being rather complex, we focus on a summation of the previous relation on legs q. This 
yields a recursion relation on the arm statistics, in other words on the choice statistics 

oo 

n D {N c ) = ^2ir D (q,N c ). (III. 11) 

9=1 
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We denote this operator by A and write it in matrix form, where rows and columns are indexed by N c , as 



A 



i 


11 


1 1 


1 1 


1 1 


1 1 




z 


2 2 


2 3 


2 4 


2 5 




§1 


1 


1 


1 


1 






2 




4 


5 









11 


_L 


_L 


1 






9 9 


o 


, | 















1 1 


1 


1 






2 3 


4 


5 















1 1 

2 4 


1 
5 


















1 1 

2 5 




\ 















(111.12) 



This operator preserves the mean number of choices (N c ) — 2. The fixed point is 



2N n 



(III. 13) 



r e 2 N °(N C - 1)!' 

The first few values of tt(N c ) are listed in table II Vl where they may be compared with Monte Carlo simulation results 
(see paper II2&). The values peak at N c = 2 and decrease rapidly for large N c . 
The mean-field entropy in eq. (|III.9(1 becomes 

1 °° 
77% 51 



MF 

— 



2N r -l 



e ^2^(^-1)! 



\ogN c ~ 0.598. 



(III. 14) 



This mean-field value is satisfactorily close to the numerical value obtained by Monte Carlo simulations (see paper 
Il2&), namely a x ~ 0.568, and is well below the upper bound log 2 = 0.693. 



TABLE IV: The first values of the limiting distribution 7v(N c ), obtained both in the mean-field approximation eq. 1111.1311 and 
numerically, by Monte Carlo simulations (from paper 





1 


2 


3 


4 


5 


tt(N c ) (Mean-field) 


0.30 


0.45 


0.19 


0.04 


0.007 


n(N c ) (Numerical) 


0.26 


0.51 


0.21 


0.03 


0.001 



3. Spatial correlations 



Spatial correlations play a crucial role in two-dimensional finite codimension tilings since they are critical systems 
with long-range correlations 21 . In particular, they condition the diffusive behavior of de Bruijn lines, they constrain 
the successive steps in the iterative construction process, and they reduce the entropy slightly below the mean-field 
value. We can calculate both short- and long-range correlations among p and q within our mean-field theory. We also 
note some expected properties of perpendicular space height correlations. 



1. Short-ranged correlations 



Consider the short-ranged correlation between the numbers of arms and legs p and q on individual vertices. Within 
the mean-field approximation, {pq)o = J2 P q 7r D(q,p)pq- If, at step D, a vertex v has p arms and q legs (see Fig. |SJ, 
then its first son vi has on average (2 + 3 + . . . (p + l))/p = {p + 3)/2 arms independently of its number of legs, and 
(1 + . . . + q)/q = (q + l)/2 legs independently of its number of arms. Similarly, its second son v-i has on average 
(p + l)/2 arms and (q + 3)/2 legs. Therefore 

(pq) D+1 = ^ D{qi p)\- — — + - — — \ (HI.15) 

p,q x ' 

= i(( M ) D + ll), 
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since (jj)d = (q)d = {N c )d — 2. The limiting value is (pq) oo = 11/3. Moreover, distribution i|III.13|) leads to 
= ((p 2 )oo — (p)D 1 ^ 2 = \/3/2 and Ap = Aq. Therefore the mean-field covariance of p's and g's is 

co V ( P) g ) = (w)o0 Ap ( f g Cc(<Z)o ° = -4/9 ~ -0.44, (III. 16) 

also in good agreement with the numerical value cov(p, q) ~ —0.36 (see paper The main contribution to those 

correlations comes from the asymmetry caused by the edge between v\ and «2- 

We can use the covariance to improve on our original "upper bound" entropy estimate eq. I|III.7|I of log 2 = 0.693. 
Whereas the original estimate arose from the mean number of choices in single steps, we can instead consider the 
geometric mean number of ways to make two steps. The product pq counts the number of two-step paths passing 
through a given vertex. Hence an improved entropy estimate is log ^Jll/i = 0.649, which is indeed closer to the 
numerical value of 0.568, though nowhere near as close as our full mean- field value in eq. (|III.14|1 . Ideally, we would 
determine the full distribution of choices (i.e. number of arms) in two-step paths, 7r 2 (-/V^, TV,?) then we could replace 
our MFT estimate eq. (fill. 141 with 

oo 

^ /F2 = E M^,N!)log^/N^, (111.17) 

which should improve upon the value in eq. (|III.14|I just as log ^/ll/3 improved upon log 2. Unfortunately, it is 
considerably more difficult to obtain ^(iV^, N%) than n(N c ). 



2. Long-ranged correlations 

To calculate long-range correlations using our iterative construction method, two phenomena compete: dc Bruijn 
lines running through a tiling destroy vertex-to-vertex correlations, while they increase distances between vertices. 
Consider two vertices separated by a distance u n on an n — > 2 tiling. The addition of successive families of lines to 
get higher codimension tilings (D — > 2 with D > n) make the tiling "swell" homogeneously, and the distances grow 
like D. When D increases, the distance between the two vertices increases so that 

u D /D~u n /n. (111.18) 

The idea is now to track such pairs of vertices, and to estimate how their correlation functions evolve at the same 
time they are moving apart. 

To track these correlation functions as D increases, we estimate correlations between the number of choices at 
vertices, N c (v). Define the correlation function between vertices v a and v b , widely separated in space, 

C,(/) = — V .V,(r„).V,(n,')- I —> XC) I ={N c (v a )N c (v b ))r-{N c ) 2 . (III. 19) 




The first average runs over the Np pairs of vertices of n — > 2 tilings separated by a distance r. The second average 
runs over all N v vertices in the tiling. In the following, we denote (N c (v a )N c (vb)) r = -f n . Moreover the mean number 
of choices (N c ) = 2, so that 

<7 n (r)= 7n -4. (111.20) 

When a new line passes through either v a or u^, their correlation diminishes. Each vertex has a probability 2/n 
that a new line goes through it (see eq. (|II.14|l '). Therefore the probability that either one of the pair (v a , Vb) will be 
affected by a new line is 4/n. A fraction 1 — 4/n of pairs separated by r are unaffected, so their correlation remains 
equal to 7„. The remaining pairs are affected, and we denote their new correlation by 7^ which we now calculate. 

Each pair (v a , Vb) that is affected gives birth to two pairs. One of its vertices, say Vb, remains unchanged while the 
other one, v a , gives birth to two vertices, a "left-son" v a i and a "right-son" v a 2 (see Fig. |SJ). On average, 

N c (v al ) = 1 + 2 + -- + N M = Nc{Va)/2 + 1/2 (m.21) 

N c {V a ) 

and N c (v a 2) = N c (v a )/2 + 3/2, since a right-son has one more arm than a left-son. Therefore, averaging on the whole 
tiling, we find 

i n = l/2({N c (v al )N c (v b )) + (N c (v a2 )N c {v b ))) = j n /2 + 2. (111.22) 
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Now we estimate j n +i, the correlation averaged over all tiling vertices after the introduction of the lines of family 
n+ 1. For original vertices not touched by the new lines, the correlation remains "f n . For original vertices that are 
touched by the new lines, the correlation becomes 7^. Therefore, 

7n+i = (1 a 4 -4/S+ 8 8/" 7 " ra (1 - 8/n)ln + 8/71 7 » + ° {l/n) - (IIL23) 

Hence 

C n+1 (u n+1 ) = (1 - A/n)C n (u n ). (111.24) 

Now iterating from n — > D we write 

D 

C D (u D ) = C n {u n ) Y[ (1-4/fc). (111.25) 

fc=n+l 

For large n and I?, 

D-l D-l 

log C D = log C n + Iog(l - 4/fc) ~ log C„ - 4 ^ ~ log C n - 4 log(D/n) (111.26) 

k—n k—n 

from which it follows that 



Using eq. (|III.18|) we have 



and finally 



D 4 C D (u D ) ~ n 4 C„(un). (111.27) 



( MD ) 4 C D ( WD ) ~ K) 4 C„(un). (111.28) 



Coo (r) ~ const, -r. (111.29) 



This expression holds only for large values of r, since we do not take into account the discreteness of tilings at very 
short range. Correlations fall off rapidly at large distance. 

3. Height-height correlations 



In two dimensions, the random tiling hypothesis eq. (|II.lJl imples logarithmic behavior of the two-point height-height 
correlation function^ 

G{r) = ((h(r) - h{0) 2 ) = -L logr (111.30) 

ttK 

in the large r limit. We give here the generic form, without regard to the specific components of our D — 2 dimensional 
height function or the particular irreducible components of the phason strain tensor K. The only requirement is that 
the elasticity tensor be positive definite, i.e. K > 0. Presumably eq. I|IH.30|) will hold for any finite D. 

We can easily derive an alternative formula for G(r) in the limit r -C D. Consider the shortest path joing two 
tiling vertices, and let s be the number of tile edges on this path. Provided s <C D, a typical path contains s distinct 
edge orientations ±e| and its length r is proportional to s. In the perpendicular space, h(s) — h(0) contains s distinct 
terms equaling ie 4- , the remaining D — 2 — s components equaling 0. Since the products (|A.6J) are of order at most 
1/D whereas (e 4 -) 2 = 1, (h(s) — h(0)) 2 ~ s when s -C D. Consequently 

G(r) oc r (111.31) 

for r <C D. Now taking the D — > 00 limit, and comparing eq. I|III.31(I with the large r limit (|III.30J| . we see that K 
must vanish in the limit of large D or else G{r) will lack monotonicity. Section llV 21 provides further evidence for the 
vanishing of K . 
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4. Finite D corrections 



The mean-field calculations of section IIII 21 neglect the fact that at the step D of the iterative process, only a 
diminishing fraction of of the vertices are visited by paths. Correcting this oversight does not alter the fixed point of 
the process, but predicts finite D corrections to it. Indeed, at step D, a fraction 2/D of the vertices are visited. Each 
of them has 2 sons and equation (|III.10|) applies only to a fraction of new vertices of order A/D when one iterates at 
step D. The remaining vertices are unaffected. Hence the correct evolution operator is 

A D = ^A + (l - j^J Id. (111.32) 

where A is the operator introduced in eq. (|III.12J) . For any finite D, the fixed point of A remains the fixed point of 
Ad- Any eigenvector e of A associated with the eigenvalue fi evolves as follows: 

A D A D - 1 ...J a .e=f[(^+0.-h)e =n(l + fc^)e ^^^e, (111.33) 

k=2 V J k=2 V 7 

where C is some constant of order 1, using the same calculation scheme as in equations (|III.26(I and (|III.27(1 . 

Hence we get power-law convergence towards the fixed point instead of exponential convergence. The spectrum of 
A can be calculated numerically, by taking finite but larger and larger approximations of A. The spectrum converges 
rapidly towards {1, 1/2, 1/3, ... , 1/p, . . .}. Eigenvalue [i = 1 is the desired fixed point. The slowest mode corresponds 
to /i = 1/2 and to a 1/D 2 behavior of the finite D corrections to ir(N c ), and therefore to the entropy (|III.9|I . Thus 
finite D corrections to the entropy associated with path-counting are expected to behave like 1 / D 2 in the mean-field 
approximation. 



IV. DE BRUIJN LINE INTERACTIONS 



The mean-field theory addresses the statistics of isolated de Bruijn lines influenced by the statistics of the D — 1 — > 2 
tiling they meander on, but with no other interactions. A crucial fact about the de Bruijn line construction of tilings 
is that lines within the same family may never cross, resulting in an effective repulsive interaction at contacts^. We 
now analyze the effect of interactions on finite size corrections and on phason elasticity. 



1. Entropy loss due to contacts in strain- free tilings 



We call "contact" between two neighbor lines of a same family a configuration where these lines are adjacent. Such 
a situation is represented in figure 01 where we have circled one contact between neighboring paths (or de Bruijn 
lines). If we follow the lines in the direction t, the local initial angles 9 and 8' they make with their global direction 
8 after the contact are constrained by the noncrossing relation 6 < 6'. This constraint divides the number of allowed 
configurations by a factor of order 2 as compared to the free path case, and reduces the total entropy by about log 2. 

To estimate the global entropy loss due to contacts in an infinite tiling, we must calculate the density of contacts. A 
crucial quantity is the distance Li between consecutive contacts with neighboring lines of family i. Inspecting Fig. 0] 
we see that contacts between lines of family i = D result from diffusion of the lines in the direction e D . Note that Li 
depends on both Zj (because wider separation reduces the frequency of contacts of diffusing lines) and fa (because the 
diffusion constant has corrections in the fa; these corrections are quadratic because of the symmetry of the problem 
through any line). In strain- free tilings, which have maximal rotational symmetry, Z, = I* — D/2 and fa = and 
Li = L is independent of i. 

Define the diffusion constant T> so that the lateral displacement u after traveling a distance L in the direction 5, 
is u 2 ~ T>L. We presume that V takes a finite limiting value in the large D limit. Thus VL ~ (l*) 2 , and there is 
approximately one contact among lines of family i in every region of area A = l*L = (Z*) 3 /T>. Adding up the total 
number of contacts expected among lines of all orientations, we estimate the loss of entropy per unit area as 

Aa ra ~w log2 = _8Plog v- (IV - 1} 

This quantity should be suitably rescaled if one is interested in the entropy per tile. The strain-free loss of entropy 
falls off like 1/D 2 because as D grows large, lines of like orientation rarely contact each other. 
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2. Phason elastic constant tensor 

We now investigate the phason strain dependence of the entropy and test the generic form assumed in cq. I|II.1|) . A 
thorough discussion of phason elasticity requires a group-theoretical analysis of rotational invariants^, each leading 
to an independent phason elastic constant, which lies beyond the scope of this paper. Rather, we will assume a generic 
form of strain and estimate its impact on the entropy. In particular, we shall see that K decays like 1/D when D 
becomes large. 

To start, consider expanding <td(_E) to second order in fa and Sli = U — I*. Note that this expansion does not 
couple explicitly angles fa and line spacings Sli because a symmetry through any line changes fa into —fa whereas 
it does not affect Sli. Therefore this expansion is the sum of a quadratic form in the <^'s and of a quadratic form 
in the Sli's. Some phason elastic constants are related to the angles fa only. The remaining ones are related to the 
variations of the line spacings k only. 

At present we consider only the k strain. To estimate the strain-dependence of the loss of entropy by contact we 
must calculate the density of contacts as a function of strain. We impose a global phason strain by varying only the 
mean distance ^ between de Bruijn lines of family i. Let 

li = l*(l + \6i) (IV.2) 

where A is a small parameter that controls the strength of the phason strain. Here Si is an arbitrary but bounded 
function of i whose mean and standard deviation obey (S) = and (S 2 ) — 1, where the averages are taken over families 
i. Owing to the bounded tile fraction constraint eq. I|II.21J| . whatever the strength of the phason strain, k is of order 
D and A is of order at most 1. The requirement of vanishing mean is imposed by the normalization condition (|II.18(I . 

Consider a pair of adjacent lines of family i, and evaluate Li, the mean distance between contacts of these lines. Thus 
T>Li ~ If, and there is approximately one contact among lines of family i in every region of area Ai = hLi = If /D. 
The entropy loss <|IV.1|) becomes 

Act w -VD(r 3 } log 2. (IV.3) 

Now consider the average (/~ 3 ). For vanishing phason strain (A — > 0), eq. (|II.15|I reveals that U = I* = D/2. 
Expanding l~ 3 in eq. (|IV.2|I for small A, using the known mean and standard deviation of Si, we obtain 

(r 3 )« (|) 3 (1 + 6A 2 ). (IV.4) 

Substituting this expression into eq. HIV.3|) yields 

Act « -(-5 + -^) 8© log 2. (IV.5) 

We find two terms, each falling off as 1/D 2 . The first term, which arises in the absence of strain, reproduces eq. (|1V.1|I . 
In a strained tiling with A ^ 0, owing to eq. I|IL21(1 . ?j remains of order D, and the loss of entropy is also of order 
1/D 2 . Although the frequency of contacts increases due to the strain, it still falls off for large D. 

Because the entropy is independent of the phason strain at the large D limit, we again conclude that the phason 
clastic constant vanishes. We have seen (eq. H11.230 ) that for bounded tile fractions the phason strain E 2 decreases 
like 1/D as D increases. Comparing eq. (|II.1(I with the strain-corrections of order 1/D 2 that we can calculate, we 
conclude that the phason elastic tensor K falls off as 1/D. 

We can understand the vanishing phason strain correction by the following qualitative argument. Consider a domain 
of the tiling of diameter A, large as compared to the tile size 1 but small as compared to the mean distance between 
two lines of a same family, 1 « A « ft For example, one can choose A = \[D. In such a domain there is at most 
one de Bruijn line per family, and no lines at all of most families. Since A tends to infinity with D, the area of the 
domain tends to infinity and the number of lines and tiles in the domain also tend to infinity. 

Locally, as D and A go to infinity, the tiling looks like a large codimension D' — > 2 tiling (where D' ~ A), with 
only one line per family (see Fig.^J). The entropy per tile in this domain tends towards a quantity o~\ oc which depends 
neither on the strain nor on the domain under consideration and is directly related to a simplified problem with one 
line in each family. The global entropy per tile Coo is the average over the whole tilingii of the local entropy density 
o~ioc- Therefore Uoo = o~i oc is independent of the de Bruijn line separations U, of the angles fa and of the global phason 
strain. 
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FIG. 9: A domain of diameter A which contains at most one de Bruijn line per family. 



V. TILINGS OF DIMENSION OTHER THAN 2 



This section discusses briefly cases of dimension other than d = 2. The high codimension d — 1 case is especially 
interesting because it is exactly solvable. It serves as an excellent illustration of several of the concepts in the paper 
and contrasts interestingly with the high codimension d = 2 case. We then present results and conjectures for the 
high codimension d = 3 case and the codimension 1 case with d = D — 1 for arbitrary d. 



1. D — > 1 tilings 

One-dimensional high codimension D — > 1 tilings are exactly solvable, as we show it in this subsection. The special 
case D = 2 was previously discussed, for example in 5 . These tilings comprise D different species of tiles randomly 
arranged on a line. This case is unusual since the entropy per tile tends to infinity with D. Still, this example is quite 
instructive since it already contains a characteristic feature of large codimension tilings: when D goes to infinity, we 
prove that the entropy no longer depends on the relative tile fractions provided these quantities remain bounded. 

The number of tilings is given by the multinomial coefficient: 

Introducing tile fractions Xi = ki/{k\ + . . . + fco), we get the entropy per tile at large D 

D 

a - 1 =-J2 Xi log Xi . (V.2) 

i=l 

This entropy is maximal at Xi — 1/D where it equals \ogD. In one dimension, for the sake of convenience in the 
calculation, we set the upper bound Xi < a/D where a > to prevent some tile species from being dominant. Then 
— \ogXi > \ogD — \oga, and hence — Yli=i x i ^°& x i — l°g-D — logo because Yli=i x i = 1- Therefore 

log.D- logo < a ^ 1 < log D. (V.3) 

The entropy grows asymptotically like logD + 0(1), independently of the exact tile fractions provided they remain 
bounded. 

Large D phason elastic constants can also be calculated. Near the entropy maximum at xi = 1/D, 

dxf 



Ki = ^a= -D. (V.4) 



Let Xi = -jy + Sxi , so that a = <r max ~ %^2, • ^ ne average gradient of the function h (representing the one- 
dimensional membrane in the D-dimensional space) is characterized by the average angle ip this membrane makes 
with the reference direction (1,1,..., 1). Specifically, (Wh) 2 = tan 2 ip ~ ip 2 . Now, 

Dj> 2 d-f$> 2 - (V.5) 

i=\ j i=l 
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Finally, (Vh) 2 = D^Sx 2 and 

a D ^^a max -^-(Vh) 2 , (V.6) 

where the phason elastic constant K = 1 is vanishingly small as compared to the entropy cr ma]t . Note the value of K 
is independent of D, and of course it agrees with the previous valued derived for D = 2. 



2. D -> 3 tilings 

Unlike the cases of d = 1 and d = 2, we have neither exact, nor extensive numerical information available for high 
codimension d — 3 tilings. Table HT1 collects available data on d = 3 strain-free boundary tilings. For D < 3 the 
entropy vanishes. No data is available for D = 5. For D = A and 6, data is available from computer simulations^^. 
The 6 — > 3 case^ deals with rhombohedron tilings with icosahedral symmetry. 

To estimate the large D limit, we insert oriented membranes into d — 3 tilings, as previously we inserted lines into 
d = 2 tilings. Membranes of like orientation are forbidden to cross. In the large D limit, we can neglect contact 
of like-oriented membranes and consider the statistics of a single membrane in generic high symmetry tilings. Let 
Md(t) denote the number of such membranes that can be inserted into tiling r and Mu be the arithmetic mean of 
this number. We find that the number of D — > 3 tilings is 

D-l 

B D = [] M D >. (V.7) 

D'=3 

analogously to eq. I|III.5(I . The entropy per tile is the logarithm of this number divided by the number of tiles. 

We obtain an upper bound on the large D entropy of d = 3 tilings in a manner similar to our upper bound of 
log 2 for for d = 2. Recall 2 D is an upper bound on the number of directed paths in d = 2 tilings, and is actually 
realized only on special tilings in which the degree Z = 4 at all sites, with two incoming and two outgoing vertices. 
Similarly, the number of directed membranes may be maximal on a special tiling in which the coordination number and 
choice numbers are most homogeneous. This is the simple cubic tiling in which each vertex is six-fold coordinated. 
The membrane orientation should be in the plane perpendicular to the (111) axis. At each vertex there are 18 
possible directed membrane configurations. This leads to Mjj < 18^ membranes containing N rhombi, a substantial 
overestimate because many of these configurations force unique membrane structures nearby. 

A true upper bound is 

M D < e^ 2N , a D ^ 3 < a 3 ^ 2 (V.8) 

because the problem of the directed membrane on the simple cubic lattice is just the 3^2 tiling problem. This value 
is an upper bound because it neglects contacts, and because generic D — *■ 3 tilings will be less homogeneous than the 
simple cubic lattice and therefore the entropy lower. 



3. D -» D - 1 tilings 

In this section, we study the entropy of strain-free D — > D — 1 tilings where the D different tiles appear with the 
same fraction 1/D and we derive rigorous lower and upper bounds: 

^<^-.<2ff=! (v . 9) 

Note that the two prefactors log 2 ~ 0.693 and 3ct 3 "^ 2 ~ 0.969 are satisfactorily close and mutually consistent. 

To derive the lower bound, we first notice that the elementary tiles of figure |21 namely hexagons and rhombic 
dodecahedrons, perfectly tile the entire two- and three-dimensional space, respectively. In fact, it is possible to tile 
any d-dimensional Euclidean space with the shadow P of a unit cube of dimension D = d + 1 projected along the 
(1, 1, . . . , 1) direction. To understand this, consider all the unit cubes of a -D-dimensional hypercubic lattice the centers 
of which belong to the diagonal hyperplane of equation x i — 0- Their projection along the (1,1,..., 1) direction is 
the desired tiling with the tile P. Now, such a polytope P can be tiled in two different ways with the D tiles of the 
D — > D — 1 tilings under consideration, as shown in figure El and independently of its neighbors. If we have a total 
of N tiles, we have N/D such polytopes, which provides 2 N / D different tilings. Thus the total number of tilings is 
bounded below by 2 N / D , which leads to the lower bound (|V.9|) on cr D ^ D_1 . 
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To understand the upper bound, we begin with the 4^3 case. In 3-dimensional tilings, we extend the notion of de 
Bruijn lines and define de Bruijn surfaces. The latter are connected sets of rhombohedral tiles which share a common 
face orientation. In 4 — > 3 tilings, there are 4 families Fi of such surfaces, since there are 4 possible edge orientations 
ef . In (i-dimensional tilings, there are (d— l)-dimensional surfaces defined similarly. Like in the two-dimensional case, 
two adjacent surfaces of a same family cannot cross but can have contacts. 

Like in the previous subsection, such a surface in a 4 — > 3 tiling can be put in one-to-one correspondence with a 
3^2 tiling, just as in a 3 — > 2 tiling, a de Bruijn line can be seen as a directed random walk on a square grid like 
a 2 — > 1 tiling. Therefore if we do not take contacts between surfaces into account, a 4 — > 3 tiling is a collection of 
independent 3 — > 2 tilings, corresponding to all the de Bruijn surfaces of a given family. Of course, this approach 
will only provide an upper bound since contacts constraint the surfaces and reduce the entropy. Now the entropy of 
3^2 tilings is known, and a fraction 3/4 of the tiles of the original 4-^3 tilings belong to the de Bruijn surfaces of 
the family under consideration. Therefore 

3 

tj < ^cr 3 " 2 = 0.242. (V.10) 

Now, we can iterate the same argument in any dimension, thus getting that 

a D^D-i D-± D _^ D _ 2 D-2 D _ 2 ^ D _ 3 3_ 

D D D v ' 

which leads to the upper bound eq. (|V.9|) . The actual 4 — > 3 entropy is strictly smaller than our bound (|V.10(1 . 
and its knowledge provides better upper bounds for larger D tilings. The best numerical estimate of this entropy is 
cr 4 ^3 = 0.214 ± 0.002 p.], which leads to a D ^ D - x < 4cr 4 ^ 3 / D = (0.856 ± 0.008)/L> for D > 4, a noticeably better 
upper bound than ljV.9l) . Similarly, the exact knowledge of any D — ► D — 1 entropy would provide a better upper 
bound for the entropy of tilings of dimension larger than D. 



VI. CONCLUSION 



This paper studied random tiling models in the limit of high rotational symmetry. We developed a mean-field 
theory for tilings based on the iterative construction of rhombus tilings introduced by de Bruijn^. Our main goal is 
to improve on the knowledge of the tiling entropy. Table lllll summarizes the progress to-datc. Evidently the mean-field 
theory is a considerable improvement upon earlier bounds. 

The relevant quantity in our mean-field theory is the number of choices at each vertex, which is related to its 
coordination number (degree). In contrast to prior work that focused on the mean number of choices, the current 
study yields the distribution of the number of choices. Table llVl demonstrates excellent agreement with numerical 
simulation. 

We argue that high symmetry tilings become insensitive to phason strain, unlike finite symmetry tilings. In other 
words, phason elasticity vanishes at large codimension. In addition to the entropy and the phason elasticity, we 
considered finite size corrections and correlation functions. 

Our mean-field theory does not take into account spatial correlations which should modify slightly the entropy. For 
example, if a vertex has a low degree, it has less chance to be visited in the iterative process, which modifies the fixed 
point ljIII.13|) and the entropy. Large degree vertices are more likely to be visited and to be split into two vertices 
of lower degree. We expect them to be less numerous than in our calculation (see table IIV|) . which accounts for the 
actual lower value of the entropy. 

Our mean-field theory could be enhanced by inclusion of correlations. Instead of single vertices (with their arms 
and legs), one might instead study the statistics of pairs of vertices or even larger tiling patches and determine how 
these statistics evolve via the iterative process, when they are split into two pieces by a path going through them. It 
would also be useful to predict how the spatial correlations calculated in the paper perturb the mean-field entropy. 

A second paper— follows this one and presents the numerical Monte Carlo study of large codimension tilings. In 
particular, it is demonstrated in this paper that in such tilings, the thermodynamic limit is restored: fixed boundaries 
do not result in the lowering of the entropy like in finite codimension cases. 
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APPENDIX A: CONSTRUCTION OF THE HEIGHT FUNCTION 



This appendix constructs^ the parallel and perpendicular spaces [£ " = R 2 and S 1 - = R D 2 ) for D — > 2 tilings as 
well as the mapping h : £ " — > £ x . We embed these spaces in the Z?-dimensional hyperspace R . Let the unit vectors 
{e^} be an orthonormal basis for R D . Their normalized projections into £" and £ ^ will be denoted, respectively, by 

ef and e^. Our first task is to represent these projections in the hyperspace basis. 
We define the two-dimensional "parallel" space fH by its orthonormal basis 

*''=\/^i> c °s(^) (A.l) 

3=0 

Given the basis for £ we define the projection operator P|| by P\\(u) = (x" • + (y" ■ and its complement 
P± = 1 — Ph. The parallel space vectors ef are the projections of suitably scaled hyperspace basis vectors, 

ef^Pu^). (A.2) 

Note that they form a regular fan 

ef =x ll cos(^) + y ll sin(^). (A.3) 

By a "regular fan" we mean a collection of equally spaced vectors arranged in the upper half pl ane. The union of 
these vectors and their negatives forms a 2£>-fold symmetric star. We set the scale factor sii = \/D/2 to normalize 

ef in eq. (|A.3j> . 

The (D — 2)-dimensional "perpendicular" space £^ is the complement in R D of £ " . Consider the projection of the 
scaled hyperspace vector s±ei into £ , which we can express as 

e+ = P±{s ± e l ) = s±(ei - P^). (A.4) 

Its magnitude can be found from 

\\Pxs±&if = 4(1 + ||P||e 4 || 2 - 2e 4 • P||6i) = 1. (A.5) 



We evaluate ||P||ei|| 2 = 1/sjj = 2/D from eq. (|A.2|> . and we evaluate e, • f^e* = \/2/Dei ■ ef = 2/D from eq. I|A.3|) 

together with eq. (|A. lfl . Then eq. I|A.5(1 yields s± = 1/ y/l — 2/D. The following result will also be useful: if i ^ j 
then 

e^ = --^ef -e]. (A.6) 

Indeed, e, ; • e 3 = S tJ = ef • ej/sjj + • e+/sj_. 

Every vertex v in a rhombus tiling occupies a position of the form 

D 

r,|(«)= 2 ef" (A.7) 

where the coefficients 7?.^ are integers. To find the values of IZi for vertex v, choose any vertex as the origin for 
then follow a path of rhombus edges from the origin to v. Adding to (or subtracting from) IZi for each edge ef along 
the path yields the values of lZi(v). The perpendicular space position of that vertex is defined as 

D 

h(r ]] ) = h(v)=Y,K i (v)ete£ ± . (A.8) 
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Values of h(r) for r in the interior of a tile can be denned using linear interpolation on the values just defined for tile 
vertices. 

The equation of a flat (strain-free) tiling, or of a locally flat tiling after coarse-graining is defined by D vectors 
rrii e £ll such that 

Tli{v) = mi ■ r||, (A. 9) 

the above equality meaning that TZi(v) is the closest integer to the right-hand side quantity. Then the (D — 2) x 2 
phason gradient tensor E = V\\h(r\\) is 

D 

E = J2^ ( A - 10 ) 

1=1 

where ej~ ® mi is a the tensor product of the column vector ef- and the line vector m^. 
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